1407.4804v1 [astro-ph.SR] 17 Jul 2014 


e 
e 


arX1v 





Astronomy & Astrophysics manuscript no. aa 
June 30, 2021 








LETTER TO THE EDITOR 


The Discovery of a Very Massive Star in W49* 


S.-W. Wul**, A. Bik!-?, Th. Henning!, A. Pasquali?, W. Brandner!, and A. Stolte4 


' Max-Planck-Institut fiir Astronomie, Kénigstuhl 17, 69117 Heidelberg, Germany 

e-mail: shiwei@mpia.de 

Department of Astronomy, Stockholm University, AlbaNova University Centre, SE-106 91 Stockholm, Sweden 
Universitat Heidelberg, Zentrum fiir Astronomie, Astronomisches Recheninstitut, Monchhofstrasse 12-14, 69120 
Heidelberg, Germany 

Argelander Institut fiir Astronomie, Auf dem Htigel 71, 53121 Bonn, Germany 


Received: 7 May 2014; accepted: 16 July 2014 
ABSTRACT 


Context. Very massive stars (M > 100 Mo) are very rare objects, but have a strong influence on their environment. The 
formation of this kind of objects is of prime importance in star formation, but observationally still poorly constrained. 
Aims. We report on the identification of a very massive star in the central cluster of the star-forming region W49. 
Methods. We investigate near-infrared K-band spectroscopic observations of W49 from VLT/ISAAC together with 
JHK images obtained with NT'T/SOFI and LBT/LUCI. We derive a spectral type of W49nrl1, the brightest star in 
the dense core of the central cluster of W49. 

Results. On the basis of its K-band spectrum, W49nrl1 is classified as an O2-3.5If* star with a K-band absolute 
magnitude of -6.27+0.10 mag. The effective temperature and bolometric correction are estimated from stars of similar 
spectral type. After comparison to the Geneva evolutionary models, we find an initial mass between 100 Mo and 180 
Mo. Varying the extinction law results in a larger initial mass range of 90 - 250 Mo. 
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1. Introduction 


Even though very massive stars (M > 100 Mo) are very 
rare, they have a strong influence on their environment via 
powerful winds and ionizing radiation, injecting large quan- 
tities of momentum and energy into the surrounding inter- 
stellar medium. Their fast evolution and the steep slope of 
the initial mass function (IMF) imply that one has to study 
the most massive star-forming regions to identify them. 
The formation mechanisms of very massive stars are by 
no means fully understood (Krumholz|/2014). For a long 
time, it was put in serious doubts whether these very mas- 
sive stars could actually form at all. Observational evi- 
dence was presented suggesting an upper mass limit of 


150 Mo (Figer||2005). However, recently 
the 


(2010) claimed existence of very massive stars up to 
300 Mo, especially in and around young massive clusters, 
such as NGC 3603, the Arches cluster and R136 in the 
Large Magellanic Cloud. 

In this letter we present the discovery of a very massive 
star in one of the most luminous Galactic HII regions: W49 
(Figfl). With dozens of OB-type stars in its core, W49 is one 
of the most important Galactic sites for studying the for- 





* Based on observations made with ESO Telescopes at the La 
Silla Paranal Observatory under programme IDs 67.C-0514 and 
073.D-0837) and on data acquired using the Large Binocular 
Telescope (LBT). 

** International Max Planck Research School for Astronomy 
and Cosmic Physics at the University of Heidelberg (IMPRS- 
HD) 


mation and evolution of massive stars in the local universe 


(Alves & Homeier|2003;|Homeier & Alves|!2005). Given its 
the Mi 


location in the plane o Way and a distance of 
11.11 ee kpe (Zhang et al.||2013), W49 is optically ob- 
scured by intervening interstellar dust. This makes an op- 
tical identification and spectral classification of the stellar 
content close to impossible, leaving the near-infrared win- 


dow (primarily K-band) for spectral classification of the 
highly obscured stars. 


Here, we report on the spectroscopic identification of 
a very massive star in W49, which we, hereafter, refer to 
as W49nrl. We first present our near-infrared observations 
(imaging and spectroscopy) of W49 (Sect. 2). The spectral 
features as well as the classification of W49nr1 are described 
in Sect. 3, where we also derive its stellar parameters like 
effective temperature (Z,g), initial mass and age. Finally, 
we briefly discuss the implications of our results in Sect. 4 
and end with conclusions in Sect. 5. 


2. Observation and data reduction 


A medium-resolution (R=10,000) K-band spectrum of 
W49nrl was obtained with ISAAC mounted on Antu (UT1) 
of ESO’s Very Large Telescope (VLT), Paranal, Chile. J- 
and H-band images were obtained with SOFI at the New 
Technology Telescope (NTT), La Silla, Chile and a K- 
band image was acquired with LUCI mounted on the Large 
Binocular Telescope (LBT), Mount Graham, Arizona. 
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Fig. 1. JHK three color image of the central area of W49. 
The massive star W49nrl is indicated with a white arrow. 
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Fig. 2. Normalised K-band spectrum of W49nrl1 with the 
emission lines annotated. 


2.1. Observations 


SOFI J- and H-band imaging observations of W49 were 
performed on 2001, June 7 (PI: J. Alves) providing a 5’ x 
5! field of view with 0.288 pixel~!. The data were taken 
with a DIT (detector integration time) of 6 s and NDIT 
(number of integrations) of 5 per saved frame. The number 
of exposures for J and H-band are 20 and 15 respectively, 
which lead to a total exposure time of 600 s (J-band) and 
450 s (H-band). The spatial resolution is ~ 0.5 — 0.7”. 


The LUCI K-band data were taken on 2009, September 
29 with the N3.75 camera, providing a 5’ x 5’ field of view 
with 0.12 pixel—!. The spatial resolution of the K image 
is ~ 0.6 — 0.7”. The observations were taken with a DIT 
of 2 s and NDIT of 10. Forty-two frames were observed, 
resulting in a total exposure time of 840 s. Sky frames 
were taken at an offset positions centred at a(2000) = 
19°087"35.8°, 5(2000) = +08°50’52.7”. 

The most luminous star in the central cluster of W49, 
W49nrl (Table |1), was observed with ISAAC in the K- 
band on 2004, August 6 (PI: J. Alves), with 3 exposures 
each with a DIT of 300s. The wavelength range covered 
by the spectrum is 2.08 ym to 2.20 um. The sky frame 
and science frames were taken with an object-sky-object 
pattern, and the nodding offset between the two science 
frames was set to 20”. HR 6572, an AOV star, was used as 
the standard star to correct for the telluric features from 
the atmosphere. It was observed about one hour before the 
science frames, in the same wavelength range as the science 
observations and with an integration time of 5 s. 


2.2. Data reduction 
2.2.1. Imaging 


The J- and H-band images were reduced using the ESO 
SOFI pipeline v1.5.2. and the K data with standard IRAF 
routines (see also 
The images were dark and flat field corrected. For the 
data, a sky frame was created from the science frames by 
rejecting the brightest pixels while combining the frames 
in pixel coordinates. A sky frame for the K-band data was 
created by combining images taken at the offset position 
and rejecting the 3 lowest and 6 highest values. 

Photometry on the JH K images was performed with 
DAOPHOT (Stetson}/1987) under IRAF. Stars are detected 
with daofind with a threshold of 30 above the background. 
Aperture photometry was performed with phot in a radius 
of (1-2) x the FWHM of the PSF. For each filter a refer- 
ence PSF model was constructed by combining the PSF of 
at least 20 objects. PSF-fitting photometry was performed 
with allstar, using the PSF model to fit all objects identified 
with a 30 confidence level over the local background. 

The Kk image of W49 has severe nebular contamination 
strongly affecting the photometry of the point sources. To 
reduce the effect of the nebulosity in the K image, we first 
removed the stars by means of PSF fitting. The residual 
frame, with all the stars subtracted, was then smoothed 
with a kernel of 12 pixels, resulting in a frame containing 
only the smooth nebular emission. This smoothed frame 
is subtracted from the original frame. After that, we per- 
formed PSF photometry on the nebular subtracted image, 
resulting in a more accurate photometry. W49nrl is located 
in the center of a compact cluster and its photometry is af- 
fected by crowding from the neighbouring stars. To quantify 
the effect of the crowding, we performed aperture photom- 
etry at the position of W49nrl1 on the residuals in the psf- 
subtracted frame. This gives an error of 15.8, 6.4 and 8.6 
% for J, H and K respectively. 

Finally, we cross-matched the obtained catalogs for each 
filter to identify the sources detected in more than one 
band. We calibrated the SOFI and LUCI photometry with 
2MASS (Skrutskie et al.|[2006). The final calibration re- 
sulted in errors in the zero points of 0.0063, 0.0071 and 
0.0055 mag for J, H and K-band respectively. We did not 
find a color dependence of the derived zero points. The final 
errors of the JHK photometry, as listed in Table|1} are a 
combination of the photometry uncertainty, errors in the 
zero points and the errors due to crowding. 


2.2.2. Spectroscopy 


The ISAAC observations of W49nrl were reduced using 
standard IRAF routines. The wavelength calibration was 
performed using the Xe and Ar arc frames. After the flat 
fielding and wavelength calibration, the sky was removed 
by subtracting the frames taken at the A and B nodding 
position. The spectra were extracted using doslit and the 
different exposures are combined to one final spectrum. To 
remove the narrow Bry emission from the diffuse nebular 
emission surrounding the cluster, the background was esti- 
mated using a Legendre function, sampling a region close 
to the star, and subtracted from the spectrum. 

The spectrum of the telluric standard star was reduced 
in the same way as the spectrum of W49nrl. Before cor- 
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Table 1. Observed and derived properties of W49nrl1 


























a(J2000) (h m s) 19:10:17.43 
6(J2000) (° ’”) +9:06:20.93 
J (mag) 16.57+0.18 
H (mag) 13.47+0.12 
K (mag) 11.93+0.10 
EW(Bry) (A) 8.24 1.7 
EW(HelIl) (A) 24+ 0.7 
EW(NIII) (A) 2.3 + 1.0 
EW(N V) (A) 2.6 + 0.9 
Spectral type 02-3.51f* 
Togp (K) 40,000 — 50,000 
BC (mag) -5.2 — -4.55 
Ay (mag) 2.0" /2.6-3.5° 


100 — 180% /90 — 250° 
1.7 - 3.1x10%* / 1.2 - 4.9x 10% 


Initial mass (Mo) 
Luminosity (L/Lo) 





Notes. ‘?) With extinction law of |Indebetouw et al.| (2005). 
©) Considering other extinction laws (see text). 


recting the spectrum of W49nrl with the standard star, 
the Bry line of the standard star was removed by fitting 
a Lorentzian profile. The resulting atmospheric transmis- 
sion spectrum was used to correct the science spectrum for 
telluric absorption using the IRAF task telluric. The signal- 
to-noise ratio (SNR) of the final spectrum is ~90. 


3. Results 
3.1. Spectral classification of W49nr1 


The final, normalized K-band spectrum of W49nrl1 is shown 
in Fig. |2} The spectrum is dominated by broad emission 
lines of Bry (2.16 wm), He II (2.189 pm), N ITI (2.116 pm) 
and N V (2.10 zm). The narrow emission component of Bry 
is a residual of the nebular subtraction. The He II and N V 
lines are indicative of an early spectral type 


(Hanson et al! 
2005). The broad emission profiles imply an origin in the 


stellar wind. These properties suggest similarities with the 


spectral classes 02-3.5If*, 02-3.51f* /WN5-7 ("slash" stars) 
and WN5-7 stars (Crowther & Walborn|2011). The sum of 
the equivalent widths (EWs) of Bry and HelII can be used 


as a discriminator between these classes. For the WN5-7 
stars, the summed EWs are expected to be above 70 A, 
while 02-3.5If* stars have a total EW between 2 and 20 A, 
with the "slash" stars lying in between. The total EW of 
both lines for W49nr1 (Table [1) is (10.6 + 1.8 A), resulting 
in a classification of W49nrl as O2-3.5If*. 





3.2. Hertzsprung-Russell diagram 


Based on the classification of W49nrl1 as an O2-3.5If* star, 
we estimated To between 40,000 and 50,000 K and the 
bolometric correction (BCx) between -5.2 and -4.55 mag 
adopting the derived values for an O2f*, O3I and an O4I 
star as representative for this class 
(2011). From our HK photometry the absolute A-band 
magnitude was derived to be -6.27 + 0.1 mag by assum- 





















ing the distance of 11.11 kpc (Zhang et al./2013), applying 
the extinction law of |Indebetouw et al.| (2005) and adopt- 


ing the intrinsic color of (HW — k) = —0.1 mag for O3I stars 


from|Martins & Plez| (2006). After applying the BCx, the 
bolometric magnitude of W49nr1 was derived to be between 


-11.47 and -10.82 mag, and the corresponding bolometric 
luminosity between 1.7 x 10° and 3.1 x 10° Lo. 

We plotted the likely parameter space of W49nr1 in the 
Hertzsprung Russel diagram (HRD) as shown in Fig. |3} As 
Top and BCx are correlated, the likely location of W49nr1 
is a diagonal ellipse. The possible locations of W49nr1 in the 
HRD was estimated by calculating the luminosity for the 
three spectral types in this class (O2If*, O3I and O41), us- 
ing their TJ, and corresponding BC. To estimate the ini- 
tial mass and age of W49nrl1, the likely parameter space in 
the HRD was compared with the Geneva theoretical stellar 
evolution models [Yusof et al.|2013), 
using models with and without stellar rotation. 

From the evolutionary tracks, the initial mass of 
W49nrl was estimated to be in the range between ~110 
and ~180 Me for models without rotation and between 
~100 and ~170 Meo for models with rotation (Fig. |3} left 
panel, vertically hashed area). While the initial mass esti- 
mate for W49nrl is insensitive to rotation, the isochrones 
for the models with and without rotation for the same 
age are very different (Fig. |3| right panel). The position of 
W49nrl suggests an upper age limit of ~2 Myrs after com- 
parison with the “non-rotating” isochrones, however, con- 
sidering the models with rotation, an age between 2 and 3 
Myrs is more likely. 

As the extinction towards W49nrl is high (Ax = 2.9 
mag) the choice of the extinction law can have a large effect 
on the derived luminosity and therefore on its initial mass 
and age. To select the best fitting extinction laws, we ap- 
plied a similar analysis to the color-color diagram of W49 
as (2012) and found that the slopes of the{Cardelli| 
(1989) and laws were not 


consistent with the observations (Wu et al, in prep). 


Indebetouw et. al. (2005 
best fitting law, but also the slopes of |Fitzpatrick| 
Nishiyama et al.| (2009) and |Rieke & Lebofsky 


consistent with the observed colors. ‘Taking into account 
all the 4 extinction laws, the estimated initial mass range 


widens to 90 - 250 Ma(see Fig. [3). 





The extintion law of ) was the 






4. Discussion and future prospectives 


In this letter we report the discovery of a very massive star 
in the center of the main cluster in W49. In the following 
we discuss the uncertainties in the derivation of the stellar 
parameters and the implications for the properties of the 
central cluster in W49. We end with a suggestion for further 
characterization of this object. 


4.1. Stellar paramaters 


Our classification of W49nrl depends on the empirical rela- 
tion between the spectral type and the equivalent width of 
the emission lines as well as the calibration of K-band bolo- 
metric corrections for early-O stars based upon atmosphere 
models derived pg Sena RC) (2011). Due to the 
very few objects used in this study, it is hard to predict 
the uncertainty of this classification and a larger number of 
stars is needed to make this calibration more reliable. 

The evolution of the very massive stars is mostly gov- 
erned by their stellar wind and mass-loss rate. These input 
parameters for stellar evolution models add uncertainties to 
the estimated initial mass and age. As a comparison to the 
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Fig. 3. HRD with the possible location of W49nrl marked 
as a vertically hashed area and horizontally hashed, taking 
into account different extinction laws. The three black filled 
circles stand for an O2If*, an O3I and an O4I star respec- 
tively. Left panel: The Geneva evolution tracks (Ekstrém| 
without stellar rotation 
ine) and with rotation (dashed line) for different masses are 
over plotted. Right panel: The main sequence isochrones 


with different ages, again without stellar rotation (solid 
line) and with rotation (dashed line). 







Geneva models we use the relation between the luminosity 
and the maximum stellar mass for homogeneous hydrogen 
burners ieee aS Za} resulting in a present-day 
mass estimate of 110- 175 Mo (and 95 - 250 Mo for taking 
into account all 4 extinction laws as discussed in Sect. 3.2). 
The stellar rotation only plays an important part in esti- 
mating the age of W49nrl from the HRD, as the rotational 
models predict longer time scales for the different evolu- 
tionary phases of the massive stars. The K-band spectrum 
is fully dominated by emission lines originating in the stel- 
lar wind, hence no estimate of the rotation can be made. 
High resolution spectroscopy of possible absorption lines to 
derive its rotation is key to understand the evolutionary 
status of this extreme star as well as the cluster. By moni- 
toring the radial velocity of the emission lines, multi-epoch 
spectroscopy could probe for a possible binary nature. 


4.2. Cluster properties 


W49nrl is located in the center of the compact central clus- 
ter in W49 (Fig. [I and thus supports the theoretical ex- 
d dy 


pectation of rapid dynamical mass segregation (e.g. 
(2005 


(2009)). |Homeier & Alves ) estimate the mass 
of this cluster as 10* Mo. is suggests that W49nrl is 


located in an environment quite similar to other very mas- 
sive stars located in Comme a) It adds to the 
growing number of stars with initial masses at or above 150 
Mo, suggesting the absence of a strict upper mass limit for 
massive stars as also suggested by numerical simulations 
2011) 

Comparing the cluster mass and the derived stellar mass 
for W49nrl1 to theoretical relations between cluster mass 
and mass of the most massive star 


shows that the presence of such a massive star is consistent 
with a normal initial mass function. A full study of the 


upper end of the IMF requires a spectroscopy classification 
of many more massive stars (Wu et al, in prep). 


5. Conclusions 


In this letter we present JH K, imaging and K-band spec- 
troscopy observations of W49nrl, the brightest star in the 
central cluster of W49. According to classification criteria 
based on the equivalent widths of Bry and HelII given by 


Crowther & Walborn| (2011), W49nrl is classified as an 


O2-3.51f* star. We estimate the effective temperature to be 
between 40,000 and 50,000 K and the bolometric luminos- 
ity between 1.7 x 10° and 3.1 x 10° L@. Comparison with 
the Geneva stellar evolutionary tracks suggests an initial 
mass range of 100 - 180 Me in the case of a single star, 
relatively independent of rotational velocity. We study the 
effect of variations in the extinction law on the stellar pa- 
rameters, resulting in a large initial mass range of 90 - 250 
Mo. Estimates of the present day mass delivers similar val- 
ues. The age depends severely on rotational velocity and 
can only be constrained to less than 3 Myrs. The next step 
will be a full spectroscopic modelling of the near-infrared 
spectrum of W49nrl resulting in stricter constraints on the 
effective temperature and luminosity. Spectral modelling 
will allow us to identify possible absorption lines at other 
wavelengths, suitable for measuring the rotational velocity. 
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